Transport through waveguides with surface disorder 
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We show that the distribution of the conductance in quasi-one-dimensional systems with surface 
disorder is correctly described by the Dorokhov-Mello-Pereyra-Kumar equation if one includes direct 
processes in the scattering matrix S through Poisson's kernel. Although our formulation is valid for 
any arbitrary number of channels, we present explicit calculations in the one channel case. Ours 
result is compared with solutions of the Schrodinger equation for waveguides with surface disorder 
calculated numerically using the _R-matrix method. 
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The prediction of statistical properties of transport 
through disordered systems is one of the fundamental 
problems in mesoscopic physics 0, Q . The transmission 
coefficient, the dimensionless conductance if we are con- 
cerned with electronic devices, or in general any scatter- 
ing quantity varies from sample to sample due to fluc- 
tuations on the microscopic configuration of disorder. 
Then, its distribution over an ensemble of systems macro- 
scopically equivalent, but microscopically different, is of 
most interest Transport through quasi-one dimen- 
sional waveguides (quantum wires in the electronic case) 
with bulk disorder has been tackled theoretically yielding 
to a Fokker-Planck equation, known as Dorokhov-Mello- 
Pereyra-Kumar (DMPK) equation [3|, |j] . This equation 
gives the evolution of the probability density distribu- 
tion of the transfer (scattering) matrix M (S) when the 
length L of the system increases. The parameter appear- 
ing as time scale in the DMPK equation is the diffusion 
time across the disordered region. This leads to a natural 
assumption, the isotropic model of uniformly distributed 
random phases for M (S) I n the jargon of nuclear 
physics, it means that there are not prompt responses 
that could arise from direct processes in the system 0. 

An analytical solution of the DMPK equation, for any 
degree of disorder, was found in the absence of time- 
reversal symmetry (/3 = 2 in Dyson's scheme Ifl), by 
mapping the problem to a free-fermion model |3> |a| • In 
the presence of time-reversal invariance (/? = 1), solu- 
tions are known only in the localized [9j and metallic 
regimes . Those solutions were recently checked with 
extensive numerical Monte Carlo simulations 

[HE3 lt 

was also found in those references that waveguides with 
surface disorder are correctly described by the solutions 
of the DMPK equation in the localized regime. Inter- 
estingly, the solutions of the DMPK equation does not 
fit with the numerical results for surface disorder in the 
crossover to- nor in- the metallic regime. 

Up to now, a complete theory that include wires or 
waveguides with rough surfaces is missing. In this letter 
we show that systems with surface disorder are correctly 
described by the DMPK equation as the direct processes, 
due to short-direct trajectories connecting both sides of 




FIG. 1: Sketch of a flat waveguide of width W supporting TV 
channels with surface disorder in a region of length L. Two 
trajectories are also sketched. The trajectory 7 passes directly 
from one to the other side while trajectory 7' passes randomly. 



the waveguides, are taken into account. We introduce 
prompt responses in a global approach as explained be- 
low. The direct processes are quantified by the ensemble 
average (S) of the iS-matrix, known in the literature as 
the optical ^-matrix [T3L IT^j . Our theoretical results are 
compared with numerical calculations based on the R- 
matrix theory |l5j to solve the Schrodinger equation for 
the one-channel case. 

The waveguide with surface disorder is sketched in 
Fig. ^ It consists in a flat waveguide of width W that 
support N open modes (or channels) with a region of 
surface disorder of length L. The scattering problem is 
studied in terms of a 2N x 2N scattering matrix S which 
has the structure 



S 



r if 

t T> 



(1) 



where r (r') and t {t') are the reflection and transmission 
matrices for incidence on the left (right) of the disor- 
dered region. The dimensionless conductance T = G/Gq, 
where Go = 2e 2 /h, is obtained from S as T — tr(tt^), ac- 
cording with Landauer's formula. 

From Fig.^ it is clear that there could be direct trans- 
mission since there are trajectories, like 7, connecting 
both sides of the waveguide without any bounce in the 
rough surface although other trajectories like 7' connect 
randomly both sides. In order to verify this statement 
we perform a numerical computation for the particular 
case N — 1. We choose W = 1; length L = 100 is 
divided into 100 pieces. To implement the surface disor- 
der the ending point of each piece is a random number 
between and a constant S that measures the strength 
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FIG. 2: Dimensionless conductance T = G/Go as a function 
of E/Eo, where Eo = h 2 n 2 /2m, of a waveguide with (a) sur- 
face disorder and direct transmission and (b) surface disorder 
but without direct transmission. The insets show the respec- 
tive waveguides. In (a) the conductance shows three regimes. 
In (b) the conductance shows a localized phase only. 



of the disorder. The ending points of each piece are con- 
nected with a spline interpolation to form a smooth disor- 
dered surface. The different realizations are done choos- 
ing different random displacements. We use a reaction 
matrix based method to solve the Schrodinger equation 
|l5j | in the scattering region with zero derivative at the 
leads. The basis states can be found using a metric for 
which scattering geometry transforms to a rectangular 
region. Finally, the scattering wavefunction is expanded 
in terms of those states. The result for the dimensionless 
conductance T, for one realization of disorder, is shown 
in Fig. |2a) as a function of the incident energy E in 
units of the transverse energy Eq = h 2 ir 2 /2m. One chan- 
nel is open for E/Eq in the range 1 to 4. Within this 
range T shows three different regions and not one as ex- 
pected in quasi-one-dimensional problems with bulk dis- 
order. For low E/Eq (from 1 to 1.6 approximately) the 
conductance is small like in a "localized" regime. As 
shown in Refs. [ill the distribution of the conduc- 
tance in the localized regime of waveguides with surface 
disorder agrees with the solution of the DMPK equation. 
A "metallic" regime is obtained for energies between 2.5 
and 4 for which the conductance is close to 1. Finally, 
there is a transition region between the "localized" and 
"metallic" regimes. The existence of direct transmission 
in the waveguides with surface disorder is clear when a 
similar waveguide is bent in such a way that the direct 
trajectories are forbidden. In Fig. Etb) the conductance 
of a bent waveguide with surface disorder is plotted. It 
shows only a localized phase. The distribution over an 
ensemble of 10 elements agrees with the DMPK solu- 
tion in the localized regime (see Fig.|3J). The appearance 
of only the localized phase in the bent waveguide yields 
strong evidence that the "metallic" region in Fig.J^a) is 
due to direct processes. 

In the phenomenological model we include the direct 
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FIG. 3: Distribution of the conductance (histogram) for an 
ensemble of 10 bent waveguides with surface disorder. It 
agrees with the solution of the DMPK equation in the lo- 
calized regime (solid line). 



processes in the solution of the DMPK equation. The 
2N x 2N scattering matrix S of the system with (S) ^ 
is related to a 2N x 2N scattering matrix So with (So) = 

through mini 



So 



(S - (S) K ) 



(i2N - (sy K s) 



-tt 



(2) 



where low is the unit matrix of dimension 2N and t n can 



be chosen as t p = t„ 
satisfies t' p t' p 



/t 



where t' is a 2N x 2N matrix which 



law - (S)k(S)^ m Here, (S) K is the 
average of S taken with Poisson's kernel, being the Ja- 
cobian of the transformation The actual measured 
direct processes are quantified by (S) which is linearly 
related to (S)k (see below). Then So is a matrix that 
describes all the scattering process as S with exception 
of the direct processes. It has a microscopic configura- 
tion of disorder that gives rise to (So) = 0. Note that 
S is reduced to So when (S)k = such that (S) = 0. 
Then the properties of the conductance in the "metal- 
lic" regime of the waveguide with surface disorder can 
be obtained from the localized regime by adding the di- 
rect processes. This means that the scattering matrices 
S and So of the waveguide in the "metallic" and local- 
ized regimes, respectively, are related through Eq. J3J). 
Since the distribution of So is given by the solution of the 
DMPK equation in the localized regime, the distribution 
of S can be obtained from the distribution of So taking 
into account the Jacobian of the transformation (J2J . 
As in Eq. QJ, Sq has the structure 



So 



r t' 
to r' 



(3) 



and both, S and So, can be 
representation; for instance 



)arameterized in a polar 
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where To is a diagonal matrix whose elements are the 
eigenvalues {r „ G [0,1]} of the matrix tot\ an d u o\ 
j = 1, 2, 3, 4, are N x N unitary matrices for (3 = 2, with 

the additional conditions = Uq and = 

(i) 

for p = 1. In the isotropic model each u is distributed 
according to the invariant measure of the unitary group, 
dfj, f^o^) ■ The probability distribution of So can be writ- 



ten as 



dPo%\S ) 



Pp({ T 0n}) 



(5) 



where sq = L/Iq with to the elastic mean free path with- 
out direct processes. Pjj_£' ({Tp n }) is the solution of the 
DMPK equation |3,la>llfl| which depends on the elements 
of To only, and 



A' 



d/ipiSo) = p f 3({T 0n }) Y[ dr 0n Y[ dfJ- ( u o^ (6) 

n=l j=l 

(see Ref. Q) is the invariant measure of So where 



Pp({ T 0n}) = Op Y[ T 0a 



N 

^n^- 2)/2 - ( 7 ) 



a<b 



A similar parameterization as Eq. Q holds for S and 
for its probability distribution we have 



Pfi({ T n\) 



(8) 



where dfip(S) has a similar expressions as Eqs. (JSJ) 
and J7J) suppressing the label "0". Now, Psg\s) in- 
clude the phases and it is obtained from the equation 
dPso\s) = dP ^(S ), which gives our main result, 



' Pf3{{T0n(S)}) 



det 



(l 2J v - (S) K (S) f K 



det 



(t 2N -S{S)\ { ) 



(2N/3+2-/3)/2 



(9) 



where V^ -1 is a normalization constant. The last term 
on the right hand side is the Jacobian of the transforma- 
tion (0), known as Poisson's kernel 0,0,0]. Here, we 
are postulating that this Jacobian is valid not only when 
So is uniformly distributed but also for Sq with uniformly 
distributed random phases. However, the average of S 
that appears in the Jacobian is not the actual value of 
(S). In fact, it is satisfied the one-to-one correspondence 
(S) = (S)k + (Sf), where (Sf) is the average of the fluc- 
tuating part of S, as can be easily seen inverting Eq. (J2J 
to write S as S = (S)k + Sf. We have verified by nu- 
merically simulating DMPK Sq's, and hence S that, for 



a given value of (S)k, (Sf) coincides with (S) — (S)k- 
Note that if (S)k — (at the same time (Sf) = 0), 
S = So and P^\{t„}) reduces to Po^fon})- Finally, 
in Eq. 10, it remains to write the To„'s in terms of S 
using Eq. J5J. As we mention before, we choose t p = tL 
once we calculate t' p numerically from the measured value 
of (S); in general it can be done by diagonalizing t' p t' p . 
Then the distribution of the conductance T can be ob- 
tained by multiple integration of Eq. . An example is 
given below. 

The case = 1 with N = 1. In what follows we 
are concerned with the = 1 symmetry only and the 
index becomes irrelevant such that we can suppress 
it everywhere. In this case, So is the 2x2 matrix [see 
Eq. Q] 



So 



-vT^o"e 2 ^ 



Toe 



i(<£o+"0o) 



(10) 



where < 0o , V'o < 27r and < To < 1 . At the level of 
So, To becomes the dimensionless conductance To, such 
that tq is replaced by To everywhere. Eq. © gives 



, fC! s dTo d(j) dipo 
MS ° ] =2WTo^^- 



(11) 



where we used that p(T ) = l/2vTo- 

As we mention above, the probability distribution of So 
is given by the solution of the DMPK equation. Although 
the solution in any phase can be used, we use the one in 
the localized regime which is well known. In this phase, 
the variable x g, d efined by To = 1 /cosh 2 xo , is Gaussian 
distributed |lllo|. 



Po S n(zo) 



1 



: exp 



— (xo 

so V 



so 



(12) 



with so = — 4 In (To). Equivalent expressions to Eqs. 11U|) 
and (|llf> are valid for S without the label "0" . Those are 
used to write To in terms of S using Eq. for a given 
(S)k- We will consider only the case (S)k = wa x , where 
<j x is one of the Pauli matrices and w a complex number. 
For t' we choose t' = yl — |w| 2 diag(e l611 , e" 102 ), where 9\ 
and 02 are arbitrary phases. From Eq. J2J we get 



T, 



(1 + H 2 )VTe"' - w*e 



Q 2ir) 



W 



1 - 2w*VTe lr i + w * 2 e 2ir > 



(13) 



where r\ = <p + ip. We note that the result is independent 
of the arbitrary phases 0\ and 82 . Finally, by direct sub- 
stitution of (S) K , Eq. (H2J and Eq. JHJ into Eq. ©, we 
obtain the result for P S0 (T,<p,ip), from which the distri- 
bution of T is obtained by integration over the variables 
4> and ip in the range to 2n. Since <fi and tp always 
appear in the combination r\ = (p + ip one integration can 
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FIG. 4: (a) Distribution of T for the one channel case and 
(3 = 1 for so = 10, and intensities of the direct processes 
w = 0.2 (dashed), 0.5 (continuous), 0.8 (long-dashed). The 
crosses are random matrix simulations (10 6 realizations). For 
clarity we present only the case w = 0.5. The agreement is 
excellent, (b) Comparison between theory (random matrix 
simulations) and numerical solutions of Schrodinger equation 
as described in the text (histogram). Here, so = 33.6 and 
w = 0.93. The agreement is excellent but the experimental 
(S) does not fit the theoretical one. 



be done easily. The result is 



Ps {T,v) = — -^=<'M> 




— I cosh 1 
so 



£o 

% 2 



(i-M 2 



2VT y/T (l-T ) 



1 - 2w*VTe ir > + w * 2 e 2ir i 



(14) 



The remaining integration over the variable rj can be done 
numerically. The results for sq = 10 and w — 0.2, 0.5 and 
0.8 are shown in Fig.QJa), as well as one random matrix 
simulation only (for clarity of the figure). The simulation 
of S is through Eq. © with So satisfying the DMPK 
equation in the localized phase. We observe an excellent 
agreement. As expected P So (T) moves from localized to 
metallic (from left to right) when \w\ take values from 
zero to one. The peak in the distribution of T is proba- 
bly reminiscent of the cut off observed in the numerical 
simulations with a higher number of channels of Ref. • 
In Fig.0Jb), we compare the theory (simulations) with 
the numerical experiment of an ensemble of six realiza- 
tions of flat waveguides with surface disorder, with the 
same parameters as those of Fig. For so = 33.6, we 
have an excellent agreement for w = 0.93. The theory 
gives (S fll ) = (S fl2 ) = 0, (S f21 ) = (1 + i) x 10- 4 , 
(S f22 ) = (-0.3 + 0.4 i) x 10~ 4 , and it is satisfied that 
(S) = 0.93cr x + (Sf). However, the theoretical value of 
(S) does not coincide with the experimental one, which is 



used for (S)k, a full matrix being more realistic. 

To conclude, we obtained the distribution of the con- 
ductance for waveguides with surface disorder. The di- 
rect transmission between both sides of the waveguide 
was taken into account in a global way. The results for 
the one channel case agree with numerical simulations, 
as well as with a numerical experiment made with R- 
matrix method. When the direct processes are increased 
the distribution moves from localized to metallic. Then, 
the direct processes could be the cause of the disagree- 
ment between the Dorokhov-Mello-Pereyra-Kumar equa- 
tion and the numerical results recently obtained. 
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a full matrix; for instance (Su) 



0.51582 + 0.06476 i. 



This is probably due to the simplicity of the model we 
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